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The determination of renormalization factors is of crucial importance. They relate the observ- 
1 ables obtained on finite, discrete lattices to their measured counterparts in the continuum in a 

suitable renormalization scheme. Therefore, they have to be computed as precisely as possible. 
A widely used approach is the nonperturbative Rome-Southampton method. It requires, however, 
a careful treatment of lattice artifacts. They are always present because simulations are done at 
lattice spacings a and momenta p with ap not necessarily small. In this paper we try to suppress 
these artifacts by subtraction of one-loop contributions in lattice perturbation theory. We compare 
results obtained from a complete one-loop subtraction with those calculated for a subtraction of 
0(a 2 ). 
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1. Introduction 

Renormalization factors relate observables computed on finite lattices to their continuum coun- 
terparts in specific renormalization schemes. Therefore, their determination should be as precise 
as possible in order to allow for a reliable comparison with experimental results. One approach 
is based on lattice perturbation theory |jl|]. However, it suffers from its intrinsic complexity, slow 
convergence and the impossibility to handle mixing with lower-dimensional operators. There- 
fore, nonperturbative methods have been developed in the last years. Among them the so-called 
Rome-Southampton method ^ (or RI-MOM scheme) is widely used because of its simple imple- 
mentation. It requires, however, gauge fixing. 

In a recent paper ^ some of us have given a comprehensive discussion and comparison of 
perturbative and nonperturbative renormalization. One of the conclusions was the possibility to 
suppress the unavoidable lattice artifacts by subtracting them perturbatively. For simple operators 
this can be done in one-loop order completely by computing the corresponding diagrams numer- 
ically. While being very effective this procedure is rather involved and not suited as a general 
method for more complex operators, especially for operators with more than one covariant deriva- 
tive. A more general approach could be based on the subtraction of one-loop terms of the order a 2 
with a being the lattice spacing. The computation of those terms has been pioneered by the Cyprus 
group |Q] and applied to various operators for different actions. 

In this paper we apply this "reduced" subtraction procedure to some exemplary operators and 
compare the results with those of the complete one-loop subtraction as given in [J3j] . We investigate 
the point operators €? s = ud, (f^ = uy^ d, = ujsj^ d and (f^ v = uo^ydiov light quarks (u, d). 
The corresponding Z factors have been measured (and chirally extrapolated) at /3 = 5.20, 5.25, 5.29 
and 5.40. We used clover improved Wilson fermions with plaquette gauge action. All results are 
computed in Landau gauge. The clover parameter csw used in the perturbative calculation is set to 
its lowest order value csw = 1- 

2. Renormalization group invariant operators 

In the RI-MOM scheme the renormalization constant Z is obtained by imposing the condition 

L tr (r R (p)r B U P )) = i (2.i) 

at p 2 = /I 2 , where T is the corresponding amputated Green function of the operator & under study. 
The Z factors relate the renormalized to the unrenormalized Green functions 

T R (p)=Z q l ZT(p), (2.2) 

with Z q being the quark field renormalization constant determined by 

7( , K{-iLxYlsHaPl) a S~ l (p)) 

Z i(p) = ■ 2/ n ' < 2 - 3 ) 

12 La sm z (ap A ) 



(S is the inverse quark propagator). Using (2. 1) we compute Z from 



^ 1 Z^tr(r(^)r Bo 1 rn (p)) = l. (2.4) 
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For operators which carry at least one space-time index (i.e. the corresponding 0(4)— multiplet 
has dimension greater than 1) we use an averaging procedure as described in [§]. 

We define the so-called renormalization group invariant (RGI) operator, which is independent 
of scale M and scheme J?', by [0] 



? rgi _ AZ ^( M ) 0*(M) = Z RG \a) ^ bare (2.5) 



with 



and 



,r,„ 1 _A«rWV <W2fc> If'*' ./r'd'] » 



AZ-(M)=(2ft^fj exp U OUpffi + ftU (2.6) 



Z RGI (a) = AZ- y '(M)Z b ^. e (M,a) . (2.7) 

g* 5 "^ 7^ and jS 1 ^ are the coupling constant, the anomalous dimensions and the /3 -function in scheme 
5?, respectively (70 and j8o are scheme independent and denote the corresponding lowest order 



coefficients). Relations ( |2.5[ ), ( |2.6| ) and ( |2.7[ ) allow us to compute the Z factor of the operator G 
in any scheme and at any scale we like, once Z RGI is known. Therefore, the knowledge of Z is 
very useful for the renormalization procedure in general. Ideally, Z RGI depends only on the lattice 
spacing a. Computed on a finite lattice, however, it suffers from lattice artifacts. For a precise 
determination it is essential to have these discretization errors under control. 

Most quantities on the lattice are computed within the so-called RI'-MOM scheme. However, 
being not covariant, this scheme is not very suitable for computing the anomalous dimensions. 



Therefore, we replace (2.7) by 



Z RG \a) = AZ^(M = ll p )z£,_ mom (M = ^)Z^; MOM (^ p ,a) . (2.8) 

For the intermediate scheme 5? we have chosen a momentum subtraction scheme. On a lattice 
with linear extent L, the scale \ip should fulfill the relation 

l/L 2 <A2 CD « jU 2 <<1/a 2 5 (29) 

then Z RGl (a) would be independent of ji p and from the resulting plateau we could read off the 
corresponding final value. The formula which is used to compute the perturbative conversion factor 
Zjy,_ M0M (p) is given in ^ together with all needed coefficients of the /3-function and anomalous 
dimensions. We will not give them here - the reader is referred to this reference. 



3. Perturbative subtraction of order a 2 

As shown in [J3j] the complete one-loop subtraction of lattice artifacts results in a very weak 
/^-dependence of the Z RGI which allows a rather precise determination. In the absence of this 
procedure there is the question whether a "reduced" subtraction can do a similar job. It could 
be based on a one-loop calculation including all possible 0(a 2 ) terms performed by the Cyprus 
group [Q]. Recently, it has been demonstrated that for some selected operators and actions the 
subtraction of those terms shows encouraging results [Q]. 
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z["_i (p,a). Then we define the subtracted Z factor as 



Let us denote the 0(a 2 ) part of the one-loop contribution to the renormalization constant by 



7 RI'-MOM 
J bare 



MC,sub 



7 Rf-MOM/ \ 2 7 {a 2 ) / „ 

'bare (P,«)MC ~ g^l-loo^KP^ 



(3-1) 



where g+ can be chosen to be either the bare lattice coupling g or the boosted coupling g B defined 
by g\\ = g 2 /P{g), P{g) is the measured plaquette at g. The final renormalization group indepen- 



dent Z factor is then computed from ( |3.1[ ) using (|2.7|), where we expect slightly different numbers 
depending on the choice of coupling g*. As suggested by the results in [^] we choose g+ = g B . 
The subtraction terms Z[_( (p,a) can be calculated to a very high precision. Therefore, the only 
significant errors to Z^ rf T MOM (j9,a)MC,sub are due to the Monte Carlo simulations. In Figure |l] 
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Figure 1: Unsubtracted and subtracted renormalization constants for the vector operator 6 (left) and the 
tensor operator T (right). 

we show the effect of subtraction (complete and 0(a 2 )) for the vector and tensor operators. The 
complete one-loop subtraction results in a clear plateau for both Z RGI factors. Using the 0{a 2 ) sub- 
traction there remains a more or less pronounced curvature. For small p 2 both subtraction methods 
agree, as they should. 



4. Fit procedure 

Compared to the complete one-loop subtraction we expect Z^ r£ T MOM (p,a)MC,sub as computed 
from (3T) to contain higher p (n > 2) terms constrained only by hypercubic symmetry. Therefore, 



we parametrize the subtracted data for each j8 in terms of hypercubic structures as follows (see 



(2-8)) 



^-MOM(^)^ MOM (^^)MC,sub = Z RG \a)/AZ y \p)+ Cl a 2 S 2 + c 2 a 2 S 4 /S 2 

+ c 3 a 2 S 6 /S 2 + c 4 a 4 S 2 + c 5 a 4 S 4 (4.1) 
+ C6 a 6 Si + c-i a 6 S 4 S 2 + eg a 6 ^6 , 

with S n = Y^i P"- The parameters c\ - cs describe the lattice artifacts. Together with the target 
parameter Z RG1 (a) we have nine parameters for this general case. In view of the limited number of 
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data points for each single /3 value (5.20, 5.25, 5.29, 5.40) we apply the ansatz (4.1) to all j8 values 
simultaneously with 

Z RG \a)/AZ y \p) -> Zf G \a)/AZf(p) , (4.2) 

where k labels the corresponding j8 value. The parameters c, are taken to be independent of j8 . This 
enhances the ratio (number of data points)/(number of fit parameters) significantly ! 

Of course, there is a certain degree of freedom in the fit procedure. One choice regards the 
interval p^ in < p 2 < p\ mx used for the fit. Inspection of data and of the results in [Q] suggests to 
use p^ nin = 10 GeV 2 for all j3 values and all considered operators. For p 2 nax we choose the cor- 
responding maximal available momentum. Another interesting point is to investigate whether the 
0{a 2 ) subtraction has been sufficient to subtract (almost) all p 1 dependence. Therefore, we perform 
two kinds of fits: one with all hypercubic structures under consideration (Zf Gl (a),c\ ,C2,C3,C4,cs, 
C6,C7,cg) and one with the a 2 dependence omitted (Zf Gl (a),C4,C5,C(,,c-i,cg,). 

Additionally, the renormalization factors are influenced by the choice for r A^g. This quantity 
enters AZ (M) in (2J5) via the corresponding coupling g y (M) (for details see [J3j] ) - We use two 
values: as standard value we take foA^ = 0.73, suggested by [0|, as a second value we choose 
r o^Ms = 0-78, close to the result given in The Sommer scale r$ is chosen as r$ = 0.5fm. 
The relation between the lattice spacing a and the inverse lattice coupling j8 is given by ro/a = 
6.050(j8 =5.20),6.603(j8 = 5.25), 6.983(j3 =5.29) and 8.285 (0 = 5.40) [§]. 



5. Results 

In Figure || we show the results for our fit procedures for the choice ToA^r = 0.73. The 
Z RGI factors for the operators V and T coincide within errors well with the complete one-loop 
subtraction calculation. The renormalization factors for S and A differ more significantly. Now, 
we compare the different fit procedures using the perturbatively subtracted 0(a 2 ) contributions 
(see (pTT|)). Generally, one can state that fitting only 0(a 4 ,a 6 ) structures still leads to reasonable 
results with smaller errors. This would mean that the subtraction of one-loop 0{a 2 ) terms takes 
into account (almost) all lattice artifacts proportional to a 2 . In addition, we do not find (at least 
for the operators considered) a remarkable difference using either a single fit for each individual jS 
data set or a combined fit for all four j3 values. 

In order to estimate the quality of the fit we compute the relative difference 

SZ(p) = (Z data (/>) -Zfit(p))/Zdata(/0, (5.1) 

where Z data (» are the data for Z^ 1 ,_ MOM {p)Z^ MOM (p,a) MC , S ub- Zfit(p) is the result of the cor- 
responding fit. In Figure || we show these differences for the operators G s and G T (the other Z 
factors for different /3 behave similarly). One recognizes that the 8Z(p) are essentially in the per 
mill range. Moreover, the figures suggest that the fit to O (a 4 , a 6 ) structures seems to be sufficient 
compared to fitting all structures in (|Q|). 

In Table | we give the results for a fit to the higher order lattice artifact terms 0(a ,a ), 
where the Z RGI are obtained from each single /3 data set. The results are given in the form 
value(err\){err2), where errl is the error of the fit parameters in the applied nonlinear fit al- 
gorithm, errl denotes the change in the results if one uses the fit range 8 GeV 2 < p 2 which is not 
totally excluded. It can serve as an indicator of systematic error. 
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Figure 2: Z RGI for the local point operators S,V,A and T at foAps = 0-73. The legends denote: "comb." 



uses the fit ansatz (4.1 )+ (4.2), "single" is based on (4.1) for each /3, "all" fits the general hypercubic structure 
and "(9(a 4 ,fl 6 )" only possible a 4 and a 6 parts. The data points for "complete sub." are obtained from the fit 
procedure discussed in [[|]; they serve as reference. 
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Figure 3: Relative errors 8Z for the scalar operator G s (left) and the tensor operator T (right) for ro A-^g = 
0.73. "all" denotes the fit to the general hypercubic structure and "<9(a 4 , a 6 )" only to possible a 4 and a 6 parts. 



Our fit procedure suggests that a "reduced" perturbative subtraction with a subsequent fit of 
0(a 4 ,a 6 ) lattice artifacts leads to reliable results. This algorithm can be used for operators with 
higher numbers of derivatives where a complete one-loop subtraction is not available. However, it 
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Op. 


roAws 


7 rgi 

^ J3=5.20 


7RGI 

^ J3=5.25 


7RGI 

P=5.29 


7RGI 

^ J3=5 .40 




0.73 


0.4592(5)(1) 


0.4542(3)(3) 


0.4523(3)(3) 


0.4503(7)(1) 




0.73 


0.7267(5)(3) 


0.7342(9)(5) 


0.7368(H)(7) 


0.7528(4)(5) 




0.73 


0.7593(4)(-ll) 


0.7654(4)(-l) 


0.7674(7)(-l) 


0.7803(3)(-4) 




0.73 


0.9038(3)(1) 


0.9148(9)(21) 


0.9154(15)(9) 


0.9387(3)(2) 




0.78 


0.4699(6)(-4) 


0.4650(6)(-2) 


0.4633(5)(-2) 


0.4579(6)(-4) 




0.78 


0.7265(5)(4) 


0.7340(9)(5) 


0.7366(H)(7) 


0.7527(5)(5) 




0.78 


0.7591(4)(-12) 


0.7652(4)(0) 


0.7672(7)(-l) 


0.7802(3)(-4) 


e 1 


0.78 


0.8910(6)(13) 


0.9042(11) (26) 


0.9085(18)(11) 


0.9332(4)(7) 



Table 1: Z RG1 for the point-like operators under consideration. The results are obtained from a fit to 0(a 4 , a 6 ) 
structures and for each single j3 data set. The fit range in momentum space is 10 GeV 2 < p 2 . The results are 
given in the form value(err\){err2) where errl denotes the error from the nonlinear fit. errl is the shift to 
value if the fit is performed for 8 GeV 2 < p 2 . The shown numbers for value correspond to the (green) full 
circles in Figure 



requires a careful investigation of the prerequisites and the parameter choices for each new opera- 
tor. Finally, one should add that the computation of the 0(g 2 a 2 ) terms in a one-loop calculation for 
those operators is also challenging. 
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